SISSA/23/2004/EP 



The Universality of Penrose Limits 
NEAR Space-Time Singularities 



M. Blau"'\ M. Borunda''''='2, M. O'Loughlin^'^, G. PapadopouW-^ 

" Institut de Physique, Universite de Neuchdtel, Rue Breguet 1 
CH-2000 Neuchdtel, Switzerland 
^ ISAS - SISSA, Via Beirut 2-4, 1-34013 Trieste, Italy 
^ INFN, Sezione di Trieste, Italy 
Department of Mathematics, King 's College London 
London WC2R 2LS, U.K. 



Abstract 



We prove that Penrose limits of metrics with arbitrary singularities of power-law type 
show a universal leading u~^-behaviour near the singularity provided that the dominant 
energy condition is satisfied and not saturated. For generic power-law singularities of 
this type the oscillator frequencies of the resulting homogeneous singular plane wave turn 
out to lie in a range which is known to allow for an analytic extension of string modes 
through the singularity. The discussion is phrased in terms of the recently obtained 
covariant characterisation of the Penrose limit; the relation with null geodesic deviation 
is explained in detail. 



^e-mail; matthias.blau@unine.ch 
^e-mail: mborunda@he.sissa.it 
^e-mail; loughlin@sissa.it 
^e-mail: gpapas@mth.kcl.ac.uk 



Contents 

1 Introduction! 



2 Penrose Limits via Geodesic DeviationI 5 



2.1 The Penrose Limit and the Null Geodesic Deviation Eauationl 6 

2.2 Penrose Limits. Geodesic Congruences and Hamilton- Jacobi Equation^ 8 



2.3 The Penrose Limits of a Static Spherically Symmetric Metrid 9 



3 Example^ 11 



3.1 Schwarzschild Plane Waves and their Homogeneous Near-Singularitv Limits! ... 11 



3.2 FRW Plane Waves and their Homogeneous Near-Singularity Limit^ 14 



4 The Universality of Penrose Limits of Power-Law Type Singularities! 18 



4.1 Szekeres-Iver Metric ] 18 

4.2 NuU Geodesies of Szckeres-Iver Metrics! 20 

4.3 Penrose Limits of Power-Law Type Singularities 21 



4.4 The Role of the Dominant Energy Condition! 23 



5 Discussion! 27 



A Hamilton-Jacobi Equations. Geodesic Congruences and Adapted Coordinated 29 



B Generalisations: Brane Metrics. Isotropic Coordinates. Null Singularities! 32 



C Curvature of Szekeres-Iver Metric^ 34 



References! 35 



1 Introduction 

The results of ^ |21 El S El have led to renewed interest in the Penrose limit construction 
El IHl • The Penrose limit associates to every space-time metric Qf^i, and choice of null 
geodesic 7 in that space-time a plane wave metric 

ds^ = g^ydx^'dx" ^2dudv + Aab{u)x''x^du^ + dx^ . (1.1) 
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Here Aab{u) is the plane wave profile matrix and the computation of the Penrose limit 
along 7(u) amounts to determining the matrix Aab{u) from the original metric g^^- 

Recently, in [S] a simple covariant characterisation and definition of the Penrose limit 
wave profile matrix Aabiu) was obtained which does not require taking any limit and 
which shows that Aab{u) directly encodes diffeomorphism invariant information about 
the original space-time metric. The geometric significance of Aab{u) (and hence of the 
Penrose limit) turns out to be that it is the standard jlUl Section 4.2] transverse null 
geodesic deviation matrix of the original metric along the null geodesic 7(n). 

The relevance of this result lies in the fact that it tells us precisely which aspects of the 
original background, namely covariant information about the rate of growth of curvature 
and geodesic deviation along a null geodesic, are detected by the Penrose limit and hence 
probed by, say, string theory in the resulting plane wave background. 

In particular, as a first step towards studying string propagation in singular (and perhaps 
time-dependent) space-time backgrounds, it is then of interest to determine the Penrose 
limits of space-time singularities in general. And here we find a pleasant surprise, namely 
a remarkably universal behaviour of Penrose limits of space-time singularities. 

It had already been found for a variety of particular brane and cosmological backgrounds 
(see e.g. [SI ^3 El El) that the exact Penrose limit is characterised by a wave profile 
of the special form 



Plane wave metrics with precisely such a profile have the scale invariance (u, v) — > 
(Ati, X~^v) and are thus homogeneous singular plane waves (HPWs) [Sl^^El- Without 
loss of generality, Aab{u) can be chosen to be diagonal and, anticipating the interpreta- 
tion of the entries of Aab as harmonic oscilator frequencies, we will parametrise Aab{u) 
as 



Moreover, in jl6l 1^1 we observed that the Penrose limit of space-time singularities of 
cosmological FRW and Schwarzschild-like metrics, i.e. the leading behaviour of the 
profile Aab{u) as one approaches the singularity, is also of the above form. This led us 
to the 

Conjecture: Penrose limits of (in some suitable sense physically reasonable) 
space-time singularities are singular homogeneous plane waves with wave 
profile Aab{u) ~ -u"^. 




(1.2) 




(1.3) 



where can be positive or negative. 
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Given the above relation between Penrose limits and geodesic deviation, this amounts 
to the conjecture that null geodesic deviation shows a universal ii~-^-behaviour near such 
space-time singularities. 

The main result of this paper is the proof of this conjecture for a large class of spheri- 
cally symmetric space-time singularities known as "singularities of power-law type" or 
Szekeres-Iyer metrics jl71ll8j . These can be spacelike, timelike or null singularities, and 
encompass practically all known spherically symmetric singular solutions of the Einstein 
equations. For technical reasons we focus on the spacelike and timelike singularities. 

It is evident that, in order to be able to say anything of substance about the behaviour 
near a singularity, some supplementary energy condition has to be imposed. This is 
what we will do and, specifically, we will prove that 

Penrose Limits of spherically symmetric spacelike or timelike singularities of 
power-law type satisfying (but not saturating) the Dominant Energy Con- 
dition (DEC) are singular homogeneous plane waves of the type (|1.3() . 

Particle and wave propagation in the background exhibit a qualitatively different 
behaviour for the 'frequency squares' bounded by 1/4 from above or below. In all the 
explicit examples that had been worked out one finds that is in the range o;^ < 1 /4. 
We will show that this is indeed also the generic behaviour: 

The resulting frequency squares lo"^ are bounded from above by 1/4 unless 
one is on the border to an extremal equation of state. 

Here by "extremal" we mean, following the terminology of , near-singularity energy- 
momentum tensors saturating the DEC, and "border" refers to a border in the Szekeres- 
Iyer phase diagrams ^7j) Figure 2 in section 4.4 of the present paper. 

The explicit proof and examples of metrics displaying a different, more singular, be- 
haviour illustrate that this universal u~'^-behaviour is not simply a consequence of di- 
mensional analysis. Rather, it is the strong form of the DEC which guarantees that 
the singularity is no worse than this and, in fact, precisely sufficiently benign (due to 
the bound on the frequencies) to allow for a consistent string propagation through the 
singularity (see e.g. ^lE])- The fact that large classes of physically reasonable metrics 
with space-time singularities give rise to such a behaviour is certainly encouraging and 
perhaps somewhat unexpected. 

In section 2 we discuss the geodesic deviation approach to Penrose Limits. In section 2.1 
we describe how to define and calculate the transverse null geodesic deviation matrix. 
We also establish the equivalence of the characterisation of Aab{u) in terms of the 
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Riemann tensor of the original metric, obtained in [S], and the description in terms of 
geodesic deviation we wih use here. Section 2.2 contains some related comments on 
null congruences and solutions to the Hamilton- Jacobi equations. In section 2.3, as an 
illustration of the geodesic deviation method, and as preparation for the calculations of 
sections 3 and 4, we determine all the Penrose limits of a static spherically symmetric 
metric. 

In section 3 we apply these results to the Schwarzschild and cosmological FRW metrics, 
obtain the general Penrose limits and discuss in some detail the emergence of the 
singular HPW behaviour in their near-singularity limits. The FRW metrics in particular, 
with the freedom in specifying their perfect-fluid equation of state, will allow us to 
anticipate some of the features that will then reappear in the general discussion of 
section 4. 

In section 4 we introduce the Szekeres-Iyer metrics (section 4.1), analyse their null 
geodesies (section 4.2) and their Penrose limits (sections 4.3). In section 4.4 we supple- 
ment this by an analysis of the DEC in these models and prove the two statements made 
above. Finally, section 5 contains various comments on applications of these results to 
string theory, open questions and future work. 

Even though somewhat outside the main line of this paper, in Appendix A we elaborate, 
following the suggestion of ^2], how to use solutions of the Hamilton- Jacobi equations 
and their corresponding null geodesic congruences to construct adapted coordinates 
which are the starting point of the more traditional approach |S1 |H] to Penrose limits. 
Appendix B summarises some minor variations of the calculations of section 2.3, and 
in Appendix C we list the non-vanishing components of the Ricci and Einstein tensors 
for the near-singularity Szekeres-Iyer metrics. 

2 Penrose Limits via Geodesic Deviation 

In |2j it was shown that the wave profile Aab{u) of the Penrose limit plane wave metric 
associated to a null geodesic 7 in a space-time with metric g^i, can be obtained directly 
from the curvature tensor of the original metric. 





Here the components refer to a parallel pseudo-orthonormal frame along 7, 



7 ~ 



2E+E- + dabE'^E'' 




with E^ = du the tangential direction. 
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This relation is reminiscent of, but should not be confused with, the well-known relation 

Raubu = -A.ab{u) (2.3) 

expressing the sole non-vanishing curvature component of the plane wave metric 

df = 2dudv + Aab{u)x''J'du^ + dx^ (2.4) 

in terms of Aab{u). Indeed, the key observation of [S] in this respect was that these 
curvature components are directly related to those of the original (pre-Penrose limit) 
metric via H2.ll) . 

Equivalently, Aah{u) can be characterised as the transverse null geodesic deviation ma- 
trix |10[ Section 4.2] of the original metric, 

= AMZ' , (2.5) 
with Z the transverse geodesic deviation vector. 

The equivalence of ((231) and the characterisation (|2.1|) of Aab{u) obtained in P! is a 
standard result in the theory of null congruences (essentially the Raychaudhuri equa- 
tion). We have found that in practice this is not only a geometrically transparent but 
frequently also a calculationally efficient way of determining the wave profile Aab{u), 
and for this reason we will explain this procedure in some detail in this section. 



2.1 The Penrose Limit and the Null Geodesic Deviation Equation 

To establish the relation between 1)2. 1(1 and 1)2. 5() . we embed the null geodesic 7 into 
some (arbitrary) null geodesic congruence. Via parallel transport one can construct a 
parallel pseudo-orthonormal frame E^, A = +,—, a, along the null geodesic congruence, 

ds^ = 2E+E- + SabE^'E^ , VuE^ = (2.6) 

such that the component E^ of the co-frame Ea is 

E+ = x''d^ , E+\^ = du , (2.7) 

i.e. the restriction of E^ to every null geodesic is the tangent vector of the null geodesic. 

Infinitesimally the congruence is characterised by the connecting vectors Z represent- 
ing the separation of corresponding points on neighbouring curves and satisfying the 
equation 

Le+Z = [E+, Z] =Ve+Z- VzE+ = . (2.8) 
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In a parallel frame, covariant derivatives along the congruence become partial deriva- 
tives, 

E^V^{Z^Ea) = Vu{Z^Ea) = {duZ^)EA , (2.9) 
and since i?+ is null one has 

g{E+,E+) = ^ {VaE+)- = . (2.10) 

Hence ()2.8() implies that {d/du)Z~ = 0, and we can set Z^ = without loss of generality. 
Then, using the geodesic equation VuE+ = 0, one finds that 

VzE+ = Z\VhE+rEa + Z\VbE+)+E+ (2.11) 

and the connecting vector equation 1)2. 8|) becomes 

— = B%Z^ (2.12) 
du 

with 

B'l = {VbE+r^E^E^V^E'^ , (2.13) 
and Z^ determined by the via 

^Z+ = Z\V,E+)+ . (2.14) 
du 

For later use we note that 1)2. 10() implies that the trace of B is 

trB^Bl = V^E>l = -^d^{^x>') , (2.15) 

explaining the ubiquity of logarithmic derivatives in the examples to be discussed below. 

It follows from (|2.12l) that the transverse components Z"" satisfy the null geodesic devi- 
ation equation 

^Z- = AMZ' . (2.16) 

where 

A% = ^Bl + B-,Bl . (2.17) 

Note that (|2.16|) is just a (time-dependent) harmonic oscillator equation with {—Aa,b{u)) 
the matrix of frequency squares. 

A routine calculation now shows that 

A% = E:E^,R\p^±^xP = , (2.18) 

with R the Riemann curvature tensor of the metric g, establishing the equivalence of 
and ((231). 



Alternatively, this can be understood in terms of the standard evolution equations for 
the expansion, shear and twist of a null geodesic congruence (see e.g. ^1 Section 4.2] or 
j2Ul Section 9.2]), which are equal to the trace, trace-free symmetric and anti-symmetric 
part of Bab, 

Bab = KK^f^Pu , Pu = 9uxx^ , (2.19) 

respectively. From this point of view, the symmetry of Aab, i.e. the vanishing of the 
antisymmetric part of i? + B^, is equivalent to the evolution equation for the twist, 
and the equivalence of (|2.1|) and (|2.5jl is the content of the evolution equation for the 
symmetric part of Bab whose trace is the Raychaudhuri equation for null geodesies. 

We see from (|2.18|) that, even though Bab depends on the properties of the null geodesic 
congruence, the particular combination of expansion, shear and twist and their deriva- 
tives appearing in Aab depends only on the components of the curvature tensor and 
the parallel frame along the original null geodesic. In particular, the geodesic deviation 
matrix Aab{u) is independent of how the null geodesic 7 is embedded into some null 
congruence. 

2.2 Penrose Limits, Geodesic Congruences and Hamilton-Jacobi Equa- 
tions 

Using the geodesic deviation approach to calculate Penrose limits, as outlined above, 
is obviously a geometrically transparent and appealing way of interpreting the Penrose 
Limit and determining Aab{u). It is somewhat less economical (economical in the sense 
of introducing the least amount of additional structure) than the characterisation ()2.1() 
of Aab{u) in terms of the Riemann tensor, which only requires a parallel frame along 
the original null geodesic and not an entire geodesic congruence. However, it may 
nevertheless be a calculationally more efficient approach if one is in a situation where one 
has a natural candidate geodesic congruence (so that one does not have to construct one 
first). In this case, the calculation of Aab{u) via geodesic deviation provides a shortcut 
to the calculation of the relevant components of the Riemann tensor. 

Both these covariant characterisations of the Penrose Limit are certainly more elegant 
than the standard systematic aproach to determining Penrose Limits (HI d |H| which not 
only relies on the existence of some special (twist-free) null geodesic congruence, but also 
requires other auxiliary constructs like Penrose coordinates (i.e. coordinates adapted to 
the congruence) and the coordinate transformation from Rosen to Brinkmann coordi- 
nates. Nevertheless, this is still frequently a useful way of performing calculations, in 
particular when combined with the systematic Hamilton-Jacobi approach to construct- 
ing adapted coordinates first pointed out in ^2]) and for this reason we provide a more 
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detailed account of this construction in Appendix A. 

In practice, therefore, the geodesic deviation approach is useful if there is a natural 
geodesic congruence. Such a null geodesic congruence can be easily constructed when- 
ever one has a solution to the Hamilton-Jacobi equation 

g'"'df,Sd^S = (2.20) 

for null geodesies. Indeed, setting 

= g'^'d.S , (2.21) 

one obviously has 

x^VpX^ = g'^'^gf^'^V^d^Sd^S = ^g^" duig'"" d^Sd^S) = , (2.22) 
so that this defines a null geodesic congruence. 

In particular, whenever the Hamilton-Jacobi equation can be separated (this includes 
all space-times investigated in this paper) the null geodesic equations become first order 
and the natural null geodesic congruence is parameterised by the integration constants 
of these first order equations. 

For a geodesic congruence defined by a solution to the Hamilton-Jacobi equation, the 
equation for Bab (|2.13jl . is 

Bab = E^E^V^d.S (2.23) 
and the equation for the trace of B, (|2.15|) . is 

B\ = V^d^S . (2.24) 

Therefore B is the covariant Hessian of the HJ function S evaluated in a parallel frame 
and the trace of B is the Laplacian of the Hamilton-Jacobi function with respect to 
the space-time metric. Since i? is a symmetric matrix, the corresponding null geodesic 
congruence is twist free. 

2.3 The Penrose Limits of a Static Spherically Symmetric Metric 

To illustrate the geodesic deviation approach to Penrose limits, and as a preparation 
for the calculations of sections 3 and 4, we now show how to quickly determine all 
the Penrose limits of a static spherically symmetric metric. We start with the metric in 
Schwarzschild-like coordinates (the extension to isotropic coordinates, brane-like metrics 
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with extended world volumes, or null metrics is straightforward and is discussed in 
Appendix B), 

dnj = de^ + sin^ 0dnl_^ . (2.25) 

Taking the Penrose limit entails first choosing a null geodesic. Because of the rotational 
symmetry in the transverse direction, without loss of generality we can choose the null 
geodesic to lie in the (t, r, 0)-plane. The symmetries reduce the geodesic equations to 
the first integrals 

i = E/fir) 

= E^/f{r)g{r)-L^/g{ry , (2.26) 

where E and L are the conserved energy and angular momentum respectively. This 
defines a natural geodesic congruence, corresponding to the Hamilton-Jacobi function 

S = -Et + 19 + R{r) (2.27) 

with 

{^^Rf = gf-'E^ - r-'gL' (2.28) 
and allows us to calculate Bab- 
We first construct the parallel frame. We have 

E+ = rdr + idt + Ode , E+\., = du , (2.29) 

and we will not need to be more specific about The transverse components are 
Ea = {El, Ea), with d = 2, . . . ,d referring to the transverse (d — l)-sphere. Since there 
is no evolution in these directions, the E^ are the obvious parallel frame components 

Ea = —^ea (2.30) 
r smt' 

with ea an orthonormal coframe for The transverse SO (d) -symmetry implies 

5ia = Aia = 

A-^{u) = A{u)5,^ . (2.31) 
Moreover, because of ()2.15|) we have 

Bii{u) = V^x^(u) -{d- l)B{u) (2.32) 
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so that we only have to calculate B22{u) = B{u), for which one finds (with, say, 62 = d^j,) 

B22 = rj/ + rje^ = du iog(r (n) sme{u)) , (2.33) 



or 



BM = 5, A log(r (n) sin Oiu)) . (2.34) 



Since 



one finds 



tvB = du log I rr'^ sin'^^^ OJf{r)g{r) ) (2.35) 



Bu{u) = dulog\^r{u)r{u)^J f{r{u))g{r{u))j . (2.36) 
Now, in general, for Bab{u) of the logarithmic derivative form 

Bab{u)=5abdu\ogKa{u) (2.37) 



one has 



and therefore 



Aab{u) = 5abKa{u)-'diKa{u) (2.38) 



A., = 6-;Arsme)-'dl{rsme) . (2.39) 



An = {rrVfg) d^{rry/fg) 

^ab ^db 

In particular, for the transverse components one has the universal result 
3 Examples 

3.1 SCHWARZSCHILD PLANE WAVES AND THEIR HOMOGENEOUS NeAR-SiNGULARITY 

Limits 

As a concrete example we will now consider the Penrose limits of the D = {d + 2)- 
dimensional Schwarzschild metric^ 

ds^ = -f{r)dt^ + firy^dr^ + r^dnj (3.1) 



^The Penrose limit for certain special (radial L = 0) null geodesies in the AdS-Schwarzschild metric 
has been discussed before e.g. in j21l I22|. The general case, using the Hamilton- Jacobi method to con- 
struct adapted coordinates, was presented in |16|. For a general discussion of limits of the Schwarzschild 
metric not depending on additional parameters like L see |23| . 
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where 

m = 1 - • (3-2) 

We are assuming that D > 4 and note that evidently only for D = 4 is m the ADM 
mass of the black hole. 

In this case we have 

r^ = E^-L^f{r)r-^ = E^-2Veff{r) , (3.3) 

where K=//(r) is the usual effective potential, with respect to which r{u) satisfies the 
Newtonian equation of motion 

^ = • (3.4) 

It follows that 

r(u) L? 



A22{u) = . . . = Am{u) - -^--r^ 

r{u) r[up 

(d+l)mL2 



r{u) 

where r{u) is the solution to the geodesic (effective potential) equation (|3.3j) . Moreover, 
since the Schwarzschild metric is a vacuum solution, this is a vacuum plane wave with 
Tr^(u) = d^^Aabiu) = 0, so that 

There are a number of facts that can be readily deduced from this result: 



• First of all, we see that the Penrose limit of the Schwarzschild metric is flat for 
radial null geodesies, L = 0. We could have anticipated this on general grounds 
because in this case the setting is SO{d + l)-invariant, implying Aab{u) ~ Sab, 
which is incompatible with TrA = unless Aab{u) = 0. This should, however, 
not be interpreted as saying that the radial Penrose limit of the Schwarzschild 
metric is Minkowski space. Rather, the space-time "ends" at the value of u at 
which r{u) = 0, say at u = 0. Perhaps the best way of thinking of this metric is 
as a time- dependent orbifold of the kind studied recently in the context of string 
cosmology (see e.g. ^25 and references therein). 

• We also learn that the Penrose limit is a symmetric plane wave (n-independent 
wave profile) if r{u) = r* is a null geodesic at constant r. Setting r = r = 0, one 
finds that 

ri-^ = {d+l)m (3.7) 
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(the familiar r = 3m photon orbit for D = 4), with the constraint 

on the ratio L/E. Precisely because they lead to symmetric plane waves, with a 
well-understood string theory quantisation, such constant r Penrose limits have 
attracted some interest in the literature. 

• Moreover we see that the resulting plane wave metric for L 7^ is singular iff the 
original null geodesic runs into the singularity, which will happen for sufficiently 
small values of L/E. 

We will now take a closer look at the u-dependence of the wave profile near the singu- 
larity r{u) = 0. We thus consider sufficiently small values oi L/E m. order to avoid the 
angular momentum barrier. 

For small values of r, the dominant term in the differential equation H3.3() for r is (unless 
L = 0, a case we already dealt with above) 

This implies that 

Thus the behaviour of the Penrose limit of the Schwarzschild metric as r ^ is 

Aii{u) = -Jis{d)u-^ (3.11) 

and 

A22{u) = ... = Am{u) = -ujls{d)u-^, (3.12) 

with frequencies 



and 

,2 2(rf + l ) 

(d + 3) 

We note the following: 



^hid) = . (3.14) 



• First of all, in this limit one finds a singular homogeneous plane wave of the type 
p. 3(1 . As we will see later, this scale invariance of the near-singularity Penrose 
limit can be attributed to the power-law scaling behaviour of the near-singularity 
Schwarzschild metric. 
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• Moreover, the dependence on L and m has dropped out. The metric thus exhibits 
a universal behaviour near the singularity which depends only on the space-time 
dimension D = but neither on the mass of the black hole nor on the angular 
momentum of the null geodesic used to approach the singularity. For example, for 
D = 4 one has 

^5(^ = 2) = ^ . (3.15) 

• The frequencies are bounded by 

u'lsid) < < ulsid) < \ . (3.16) 

• Finally, we note that the above result is also valid for (A)dS black holes since the 
presence of a cosmological constant is irrelevant close to the singularity. 



3.2 FRW Plane Waves and their Homogeneous Near-Singularity Limits 

As another example we consider the Penrose limit of the D = (n + l)-dimensional FRW 
metric 

ds" = -df + a{t)\dr^ + fk{rfdQl_^) , (3.17) 
where /fe(r) = r, sinr, sinhr for = 0, +1, — 1 respectively. 

Since the spatial slices are maximally symmetric, up to isometries there is a unique null 
geodesic and hence a unique Penrose limit. So without loss of generality we shall consider 
null geodesies which have vanishing angular momentum on the transverse sphere. 

Then, with a suitable scaling of the affine parameter, the null geodesic equations can 
be written as 

^t{u) = ±a{t{u))-^ , -^r(ii) = a(t(n))-2 (3.18) 
(and in what follows, we choose the upper sign in the first equation). Thus 

E+ = du = a-^dt + a-^dr , (3.19) 
and this can be extended to a parallel pseudo-orthonormal frame by 

E- = \{-adt + dr) 

Ea = {afkT^ea , (3.20) 
where Ca is an orthonormal frame for d^\, d = n — 1. 

The transverse rotational symmetry implies that Bab{u) = B(u)dab and Aabiu) = 
A{u)5ab- Therefore, to determine B{u) it suffices to compute the trace of Ba,b{u), 

trB = 5„log(a-Vr') ' (3-21) 



14 



implying 

B{u) = du\ogiafk) . (3.22) 

Using ^fk = -kfk, one finds 

Aab{u) = 5abA{u) = 5aA^-^^] • 3.23 

\a[u) a[up J 

This is the precise analogue of the expression 1)2. 4U() obtained in the static spherically 
symmetric case, the spatial curvature k now playing the role of the angular momentum 
L\ 

This can now be rewritten in a variety of ways to obtain insight into the properties of 
this FRW plane wave. For example, writing this in terms of t-derivatives (in order to 
make use of the Friedmann equations) , we find 

where a{t) is determined by the Einstein (Friedmann) equations, u{t) by du = a{t)dt, 
and a' = ^a. The Friedmann equations 

a'{tf + k _ WttG 



a{tY Ti{n — 1) 

a"{t) SttG 



Pit) 

[{n-2)p{t)+nP{t)] , (3.25) 



a{t) n{n — 1) 

imply 

5:fi)l±A_5|)=J=^L|,„) + P(,), , (3.26) 

a(tj^ a[t) [n — 1) 

so that one finds that the wave profile of the FRW plane wave can be written compactly 
as 



= -Z 7 TTTnI • V^-^^) 



SttG p{u) + P{u) 
n — 1 a{u) 

One immediate consequence is that the Penrose limit is flat if and only if p + P = 0, 
corresponding to having as the only matter content a cosmological constant. This is 
in agreement with the result [S] that every Penrose limit of a maximally symmetric 
space-time is flat. 

We will now study the behaviour of A{u) near a singularity, and to be specific we choose 
the usual equation of state 

P{t) = wp{t) . (3.28) 

We consider w > —1 {w = —1 would correspond to the case p + P = already dealt 
with above) and introduce the positive parameter 

h{n,w) = —^ — - (3.29) 

^ ' n{l + w) ^ ' 
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and the positive constant (constant by the continuity equation for p) 

Ch = -^^^^p{t)a{tfl^ , (3.30) 
n[n — 1) 

in terms of which the Friedmann equations read 

a'{tf = Cha{t)^'^^~'^'^/'' -k (3.31) 

a"{t) = ^C7/,a(t)(^-2)/'^ . (3.32) 

Thus the universe is decelerating for < /i < 1 and accelerating for /i > 1, the critical 
case h = 1 corresponding to Wc = —1 + 2/n (the familiar dark energy threshold Wc = 
— 1/3 for n = 3). 

We first consider the case A; = 0. In that case one has 

a{t) ~ , (3.33) 

and therefore 

a{u)r^u^/^+^ . (3.34) 
It then follows immediately from (|3.23l) that, more explicitly, the w-dependence of A{u) 

is2 

A{u) = -LO^pjiwih, k = 0)^-2 , (3.35) 

where 



^FRwi^^ k = 0)= ^ . (3.36) 



h 

[I Vhf 

We see that the Penrose limit of a spatially flat FRW universe with equation of state 
P = wp is exactly a singular homogeneous plane wave of the type (jTHI. 

The frequency square ujpj^yy{h, k = 0) has the following properties: 

• Since 

= 0) = ^Fw(l/^;^ = 0) ) (3-37) 

for every accelerating (inflating) solution of the k = Friedmann equations there 
is precisely one decelerating solution with the same Penrose limit. The self-dual 
point h = 1 corresponds to the linear time-evolution a{t) ~ t. 

• The frequency squares are again bounded by 

iolnwih,k = 0) < ^ , (3.38) 

with equality attained for h = 1. 
^This generalises the result reported in [S]. 
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• Curiously, the frequencies obtained in the Penrose hmit of the Schwarzschild met- 
ric (in ah but one of the directions) are precisely those of a dust-filled FRW 
universe, P = w = 0, of the same dimension n = d + 1, 

u;ls{d)=u;^FRw{h,k = 0) , (3.39) 

e.g. 6/25 for n = 3. 

It is clear that for A; = 0, when only the first term in 1)3. 23(1 is present, this homogeneous 
-behaviour is a consequence of the exact power-law behaviour of a{t) and hence a(n). 
Let us now consider what happens for k ^ 0, when there is a competition between the 
two terms in ()3.23|) as one approaches the singularity. 

One might like to argue that, even for k ^ 0, one finds the same behaviour provided that 
the matter term dominates over the curvature term in the Friedmann equation (|3.31() 
as a — > 0. This happens for < h < 1, and this argument is correct as one can also see 
that in this range the first term in 1)3. 23() . proportional to tt~^, indeed dominates over 
the second (curvature) term which (cf. 1)3. 34() ) is proportional to Thus for 

< h < 1 the near-singularity limit of the FRW plane wave is a homogeneous plane 
wave with /^-independent frequencies (|3.36|) . 

0<h<l: u;lRwih,k)=ujlRwih,k = 0) = j^-^j^. (3.40) 

Now let us look at what happens as one passes from a decelerating to a critical (h = 1) 
and then accelerating (/i > 1) universe. First of all, for h = 1, both terms on the right 
hand side of the Friedmann equation 1)3. 31(1 contribute equally (they are constant), 
and correspondingly both terms in (|3.23|) are proportional to u~'^. Thus one finds a 
homogeneous plane wave, but with a curvature-induced shift of the frequency, 

I 

h = l: ujj,jiw{h = l,k)=ujljiyy{h = l,k = 0) + kc'^ = - + kc^ (3.41) 

for some constant c. In particular, in the spatially closed case k = +1 (this requires 
C/i > 1), one now finds frequency squares that are larger than 1/4. This is a borderline 
behaviour in the sense that, as can easily be seen from 1)3. 31|) . the initial singularity for 
k = +1 ceases to exist for h > 1. 

It thus remains to discuss the case k = —1 and h > 1. Given the previous discussion, 
one might be tempted to think that now the second term in 1)3. 23() will dominate over 
the first, leading to a non-homogeneous and more singular u~'^^^^^~^^^-heh.SLviour. This 
is, however, not the ()3.34|) now represents the leading behaviour at large a(n). 

At small a{u), the leading behaviour is, exactly as for h = 1, determined by the constant 
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curvature term in H3.31() . Thus even in this case one finds a singular homogeneous plane 
wave, with frequency 



once again bounded from above by 1/4. 

4 The Universality of Penrose Limits of Power-Law Type Singu- 
larities 

In the previous section we have presented some evidence for a remarkable 

Conjecture: Penrose limits of physically reasonable space-time singularities 
are singular homogeneous plane waves with wave profile Aab{u) ~ 

In this section we will show how to prove this conjecture for a very large class of physical 
singularities of spherically symmetric type. We will in the process also see some examples 
of "extreme" stress-energy tensors that give rise to a different behaviour. 

4.1 Szekeres-Iyer Metrics 

The homogeneity of the Penrose limit (geodesic deviation) that we have found in the 
above examples appears to reflect a power-law scaling behaviour of the metric near the 
singularity. Thus to assess the generality of this kind of result, one needs to enquire 
about the generality of space-time singularities exhibiting such a power-law behaviour. 

In (see also |18|). in the context of investigations of the Cosmic Censorship Hypoth- 
esis, Szekeres and Iyer studied a large class of four-dimensional spherically symmetric 
metrics they dubbed "metrics with power-law type singularities". Such metrics en- 
compass practically all known singular spherically symmetric solutions of the Einstein 
equations, in particular all the FRW metrics, Lemaitre-Tolman-Bondi dust solutions, 
cosmological singularities of the Lifshitz-Khalatnikov type, as well as other types of 
metrics with null singularities. 

In "double-null form" , these metrics (in d -|- 2 dimensions) take the form 




(3.42) 



ds2 = -e^(f^'^)dC/dy + e 




where A(U,V) and B(U,V) have expansions 



A{U,V) 
B{U,V) 



p\nx(U,V) + regular terms 
q In x{U, V) + regular terms 



(4.2) 
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near the singularity surface x(U,V) = 0. 

Generically, the residual coordinate transformations U U'{U), V ^'(^) preserving 
the form of the metric (|4.H) can be used to make x{U, V) linear in U and V, 

x{U,V) = kU + lV , k,l = ±l,0 , (4.3) 

with r] = kl = 1,0,-1 corresponding to spacelike, null and timelike singularities re- 
spectively. This choice of gauge essentially fixes the coordinates uniquely, and thus the 
"critical exponents" p and q contain diffeomorphism invariant information. 

The Schwarzschild metric, for example, has p = {\ — d)/d and q = 2/d, as is readily seen 
by starting with the metric in Eddington-Finkelstein or Kruskal-Szekeres coordinates 
and transforming to the Szekeres-Iyer gauge. Alternatively, if one is just interested in 
the leading behaviour of the metric, one can simply expand the metric near r = and 
then go to tortoise-coordinates. 

We will focus on the behaviour of these geometries near the singularity at x = 0, where 
the metric is 

ds"^ = -xPdUdV + xHnl . (4.4) 

For generic situations this leading behaviour is sufficient to discuss the physics near the 
singularity. In certain special cases, for particular values of p, q or for null singulari- 
ties, this leading behaviour cancels in certain components of the Einstein tensor and 
the subleading terms in the above metric become important for a full analysis of the 
singularities jl7l I18j . The analysis then becomes more subtle and we will not discuss 
these cases here. In the following we will consider exclusively the metric ()4.4() which, 
for ?7 7^ and generic values oip and q, captures the dominant behaviour of the physics 
near the singularity. 

For 7^ we define y = kU — IV and choose k = rjl = 1. Then the metric takes the 
form 

ds'^ = Tjx^dy'^ - rjx'Pdx'^ + x'^dQ^ . (4.5) 

With the further definition r = x'^/^ (for (7 7^ 0), this has the standard form of a 
spherically symmetric metric. We will come back to this below in order to be able to 
make direct use of the analysis of section 2.3. 

For ?7 = 0, on the other hand, we choose x = U , y = —V, so that the metric is 

ds"^ = xPdxdy + x'^dnj , (4.6) 
which has the form of the spherically symmetric null metrics analysed in Appendix B. 
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4.2 Null Geodesics of Szekeres-Iyer Metrics 



In terms of the conserved momenta P and L associated with y and, say, the colatitude 
9 of the (i-sphere, 

dnl = de'^ + sin^edill_-^ , (4.7) 

in particular 

x'ie = L , (4.8) 
the null geodesic condition (for any r/) is equivalent to 

i;2 = p2a;-2p + r]L'^x-P-i , (4.9) 

To understand the null geodesics near x = 0, we begin by extracting as much information 
as possible from this equation, recalling that due to the expansion around x = we can 
only trust the leading behaviour of this equation as x — > 0. 

Unless p = q, one of the two terms on the right-hand-side of H4.9() will dominate as 
X ^ 0, and thus the generic behaviour of a null geodesic near x = is identical to that 
of a geodesic with either L = or P = 0. In the former case, one finds 

Behaviour 1: x(u) ~ n^/(^+^) (4.10) 

unless p = —1 when x{u) ~ expu. We are only interested in those geodesics which 
run into the singularity at x = at finite u. This happens only for p > —1. In the 
latter case, corresponding to null geodesics which asymptotically, as x ^ 0, behave like 
geodesics with P = 0, we evidently need r] = +1 (a spacelike singularity), which leads 
to 

Behaviour 2: x(n) ~ ^VCp+^+a) (4 

unless p + q = —2 which again leads to an exponential behaviour. These null geodesics 
run into the singularity at finite u for p + q > —2. 

For r] = +1, the situation regarding null geodesics that reach the singularity at finite u 
is summarised in the following table. 



Conditions on (P, L) 


Constraints on {p, q) 


Behaviour 


P / 0,L = 


p> -1 


1 


p = 0,L / 


p + q> -2 


2 


p/0,L/0 


p > q,p > —1 


1 


p/0,L/0 


p < q,p + q > -2 


2 


p/0,L/0 


p = q > —1 


1 = 2 



(4.12) 
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For r] = —1, the situation is largely analogous, the main difference being that now the 
second term in ()4.9|) acts as an angular momentum barrier preventing e.g. geodesies 
with L 7^ for q > p from reaching the singularity at a; = 0. These cases are indicated 
by a '— ' in the table below. For the same reason, for p = q one finds the constraint 
\P\ > \L\. 



Conditions on (P, L) 


Constraints on (p, q) 


Behaviour 


p / 0,L = 


p> -1 


1 


p = 0,L / 






p/0,L/0 


p > q,p > —1 


1 


p/0,L/0 


p<q 




1^1 > 1^1 


p = q > —1 


1 = 2 



(4.13) 



Finally, for 77 = we find Behaviour 1 for all values of p and q, with the corresponding 
constraint p > —1. 



4.3 Penrose Limits of Power-Law Type Singularities 

We will now determine the Penrose limits of the Szekeres-Iyer metrics along the null 
geodesies reaching the singularity x = at finite u. 

For r/ 7^ we notice that the metric is simply a special case of a spherically symmetric 
metric and thus can be treated using the analysis of section 2.3. Indeed, with t = y and 
r = x'^l'^ {(1^ 0), the metric (|4.5|) takes the form 

ds^ = 7]xPdy^ -T]xPdx^ + x''dn^ (4.14) 

= ^r^p/idt^ -%^(p-i+^yidr^+r^dn^ (4.15) 

where in the second line the notation of t and r is adapted to the case of rj = —1 where 
the singularity is timelike and t is time. We will continue to use this notation even for 
spacelike singularities where t is actually spacelike. 

The case g = is special, but actually corresponds to a shell crossing singularity ^7] 
which is usually not considered to be a true singularity as the transverse sphere is of 
constant radius x'^ = 1. Such singularities arise for instance for certain collisions of 
spherical dust shells. From here on we will only discuss q ^ 0. 

Referring to section 2.3 where such a spherically symmetric metric was treated, we can 
identify 

/(r) = -ryr^P/^ (4.I6) 
g(r) = -^r^(P~i+^)/i. (4.17) 
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We can now appeal to 1)2.361 12.34(1 to deduce that 



Bab = Sabdulog{Ka{u)) (4.18) 

with 

Ki(u) = r(n)r(n)2(P+i)/'? 

K2iu) = r{u)sm{e{u)) (4.19) 

As shown in Appendix B, it is straightforward to extend this analysis to the case of null 
singularities, rj = 0, with the result that the expressions for Ka{u) are identical to those 
given for rj = ±1 in ()4.19|) . We can thus treat all three cases simultaneously. 

It follows from the analysis of the previous section that the only possibility of interest 
for r{u) = x{uY/'^ is the power-law behaviour 

r{u) = u" , (4.20) 

with 

Behaviour 1: p > -1 a = q/2{p + 1) (4.21) 

Behaviour 2: p + q > -2 a = q/{p + q + 2) . (4.22) 

Clearly, then, K\{u) is also a simple power of u. Specifically one has (since we are 
interested in the logarithmic derivatives of -R'i(n), proportionality factors are irrelevant) 

Behaviour 1: Ki[u) ~ r{u) 

Behaviour 2: Ki{u) '-^ r{uf/'i . (4.23) 
Thus the corresponding component of Aab{u) is 

Behaviour 1: Aii{u) = )\ \ = a(a — l)n~^ 

Ai(n) 

Behaviour 2: Aii{u) = :^lM =palq(pa/q - l)u~'^ . (4.24) 

Ki{u) 

and the Penrose limit behaves as a singular homogeneous plane wave in this direction. 
Since b{b — 1) has a minimum —1/4 at 6 = 1/2, this leads to the bound 

^l<\- (4-25) 

This is the same range that we found empirically for both the Schwarzschild and FRW 
plane waves near the singularity. 
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The behaviour of A22 is more subtle due to the dependence of K2{u) on sin^(n). The 
general behaviour is as in (|2.4UI) . namely 

A22{u) = ^ - ^ . (4.26) 

With the power-law behaviour r{u) = u", the first term is always proportional to 
This term is dominant as n ^ when a < 1/2, while it is the second term that dominates 
for a > 1/2 (and leads to a strongly singular plane wave with profile ~ In the 

special case a = 1/2, both terms are proportional to Thus one has, for L 7^ 0, 

tul = a(l - a) < i (4.27) 

u;l = ^ + L^>^ (4.28) 

(4.29) 

When T] = 0, Behaviour 1 arises in the entire p-q plane. Thus Lof < 1/4 and a can take 
any of the three values above with the special value a = 1/2 corresponding to the line 
q = p + 1. When p > q and rj ^ 0, a = q/2{p + 1) and thus always a < 1/2. On the 
other hand, when p < q and r] = 1 we see that a = q/{p + q + 2) can take on any value, 
with a = 1/2 along the line q = p + 2 and a > 1/2 for g > p + 2. When p < q and 
rj = —1 we cannot reach the singularity along a geodesic with L ^ 0. 

When L = 0, only the first term in 1)4. 26(1 is present, and one thus finds 1)4. 27(1 for 
all values of a. Since L = implies Behaviour 1, this means a = q/2{p + 1). Along 
the special line q = 2{p + 1) one has a = 1 and one finds the "flat" Penrose limit 
Aii{u) = A22{u) = 0. In particular, this happens for radial null geodesies in the 
Schwarzschild metric {p = {1 — d)/d and q = 2/d), as already noticed in section 3.1. 

These results are summarised in Figure (la) for r] = 1 and in Figure (lb) for rj = —1. 



r{u) = a < 


1 . 

2 • 


^22 1 


[u)- 




2 


a = 


1 . 

2 • 


^22 1 


[u]- 


-> — c<j|m~ 


2 


a > 


1 . 

2 • 


^22 1 


[u]- 




-4a 



4.4 The Role of the Dominant Energy Condition 

We thus see that while we frequently obtain a singular HPW with lj'^ < 1/4 in the 
Penrose limit, other possibilities do arise. For timelike singularities, the situation is 
clear: 

Penrose Limits of timelike spherically symmetric singularities of power-law 
type are singular HPWs with frequency squares bounded from above by 1/4. 

We will now show that for spacelike singularities a different behaviour can occur only 
when the strict Dominant Energy Condition (DEC) is violated, in particular, that the 
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Figure 1: The Penrose Limit Phase Diagram in the p — g plane for (a) spacelike (77 = +1) 
and (b) timelike (rj = —1) singularities. Singular HPWs arise in the light-shaded regions 
whereas in the dark-shaded region there are Penrose limits leading to strongly singular 
(and non-homogeneous) plane waves, (a) The diagram is bounded on the left by the 
lines p = —1 and p + q = —2. The dashed line a = 1/2 q = p + 2 separates the 
two regions, and only along that line one finds singular HPWs with a;! > 1/4. (b) For 
rf = —1, one finds singular HPWs with < 1/4 for all {p,q) with p>— 1, 0^2 = 1/4 
arising only along the dashed line a = 1/2 for zero angular momentum, L = 0. 
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strongly singular region (the dark-shaded region in Figure (la)) is excluded by the 
requirement that the DEC be satisfied but not saturated. 

We begin by recalling the definition of the Dominant Energy Condition on the stress- 
energy tensor (or Einstein tensor G'^) ITU]: for every timelike vector v^, T^yV^v"^ > 0, 
and T^^v"^ is a non-spacelike vector. This may be interpreted as saying that for any 
observer the local energy density is non-negative and the energy flux causal. 

Next we recall that a stress-energy tensor is said to be of type I JOj if has one timelike 
and three (more generally, d + 1) spacelike eigenvectors. The corresponding eigenvalues 
are —p {p the energy density) and the principal pressures Pa, a = 1, . . . ,d + 1. For a 
stress-energy tensor of type I, the DEC is equivalent to 

P>\Pa\ . (4.30) 

We say that the strict DEC is satisfied if these are strict inequalities and we will see 
that the "extremal" matter content (or equation of state) for which at least one of the 
inequalities is saturated will play a special role in the following. 

The Einstein tensor of the metric ()4.5() is diagonal (Appendix C), 

id(d - l)x-'? - |r?dg((d - l)q + 2p)x-^P+'^^ 
\d{d - l)x-'' + li]dq{2p + 4 - (d + l)g)x-(P+2) 

\{d - l){d - 2)4x-« + \ri{Ap - 4g + Aqd - d{d - l)q^)S^jX-^P^'^^ (4.31) 

and hence clearly of type I. For spacelike singularities, r] = +1, we have energy density 
p = —G^, radial pressure Pr = and transverse pressures Pi = G\, while for 77 = — 1 
the roles of G% and are interchanged. 

Since for q > p + 2 the first term in G% and G^ dominates over the second term as 
X ^ 0, it is obvious that for q> p + 2 the relation between p and becomes extremal 
as X — > 0, 

Gl-Gl^{) ^ p + Pr^O . (4.32) 

Put differently, q < p + 2 is a necessary condition for the strict DEC to hold. Since 
strongly singular plane waves (the dark-shaded region in Figure (la)) arise only for 
q > p + 2, we have thus established that 

Penrose Limits of spacelike spherically symmetric singularities of power-law 
type satisfying the strict Dominant Energy Condition are singular HPWs. 

Since frequency squares exceeding 1/4 can only occur along the line q = p + 2 itself, we 
can also conclude that 



X 

Gl 
G] 
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Figure 2: The Penrose Limit + DEC Phase Diagram in the p — q plane for (a) spacehke 
{rj = +1) and (b) timehke {rj = —1) singularities. In the highlighted regions A and B the 
DEC is satisfied (but not saturated), (a) the strongly singular (and non- homogeneous) 
plane waves of the dark-shaded region with extremal equation of state are excluded, 
and sing ular HPWs with a;| > 1/4 arise only along the boundary q = p + 2 to the 
extremal equation of state, (b) the Penrose limits are sing ular HPWs with w| < 1/4, 
with = 1/4 only along the dashed line q = p + 1. 

the resulting frequency squares to'^ are bounded from above by 1/4 unless 
one is on the border to an extremal equation of state. 

A more detailed analysis of the DEC (as performed for d = 2 in il7|), shows that the 
actual region in which the strict DEC is satisfied (taking into account also the conditions 
involving the transverse pressures Pi), is more constrained. For spacelike singularities, 
this is the (infinite) region bounded by the lines 

q = 2/d, q=p + 2, g = 2(p + 1) , (4.33) 

displayed as the highlighted region A of Figure (2a) (drawn here for d = 2). A look at 
this figure confirms the results we have obtained above. 

For timelike singularities, the region where the strong DEC is satisfied is considerably 
smaller - it is a finite subset of the strip bounded by the lines q = and q = 2/d, 
indicated (for d = 2) as the highlighted region B of Figure (2b). While of no conse- 
quence for the present discussion, the fact that in region B the pressures Pr and Pi 
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cannot simultaneously be positive plays an important role in the discussion of Cosmic 
Censorship in TT. 

The fact that the Penrose limits of timelike singularities always behave as while in 
the spacelike case strongly singular Penrose limits can arise (even though only for metrics 
violating the strict DEC), might give the impression that timelike (naked) singularities 
are in some sense better behaved than spacelike (censored) singularities. We believe 
that this should rather be viewed as an indication that massless particles are inadequate 
for probing the geometry of timelike singularities since, for large regions in the {p, q)- 
parameter space, the angular momentum barrier prevents non-radial null geodesies from 
reaching (and hence probing) the singularity. From this point of view, it is much more 
significant that for spacelike singularities massless particle probes with arbitrary angular 
momentum all detect homogeneous singular plane waves provided that the strict DEC 
is satisfied. 

We conclude this section with a comment on the case of null singularities of power-law 
type (ry = 0) which are analysed in As mentioned in section 4.1, in this case 

some of the leading components of the Einstein tensor vanish and hence one (somewhat 
trivially) ends up with an extremal equation of state. Thus no interesting constraints 
arise from imposing the DEC, and using only the leading form ()4.6() of the metric cannot 
be the basis for a full analysis which is more subtle and will be left for future work. 

5 Discussion 

We have shown that space-time singularities exhibit a remarkably universal homoge- 
neous M~^-behaviour in the Penrose Limit. We have established this in complete gener- 
ality for timelike singularities of power-law type and have also shown that for spacelike 
singularities of power-law type, for which more singular Penrose limits are possible, this 
ii~-^-behaviour is implied by demanding the strict DEC. 

Perhaps the main implications of this result are for the study of string theory in singular 
and/or time-dependent backgrounds. In general, because of the simplifications brought 
about by the existence of a natural light cone gauge plane wave (and more general 
pp-wave) backgrounds provide an ideal setting for studying such problems. Now, as 
we have seen, the Penrose limits of a large class of singularities are always at least 
as singular as u^'^. Thus "weakly singular" plane waves with profile ~ n~", a < 
2, while perhaps interesting as toy-models of time-dependent backgrounds in string 
theory jl9l I26j . do not actually arise as Penrose limits of standard cosmological or 
other singularities. Moreover, a strongly singular behaviour with a > 2 can only arise 
for metrics violating the strict DEC. This singles out the singular HPWs with profile 
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~ n ^ as the backgrounds to consider in order to obtain insight into the properties of 
string theory near physically reasonable space-time singularities. 

String theory in precisely this class of backgrounds, and in precisely the frequency range 
a;^ < 1/4 that we generically find, has already been studied in some detail and shown to 
be exactly solvable in |14j . The significance of this bound on the frequencies for string 
theory (also noted in the earlier studies |19j ^ lies in the fact that the properties of Bessel 
functions allow for an analytic continuation of string modes through the singularity in 
this case (due to the standard shift by 1/4 in the Bessel equation).^ Other aspects of 
string theory in this class of backgrounds still remain (and deserve) to be explored. In 
particular, since the Penrose limit can be considered as the origin of a string expansion 
around the original background jHj (see |2H| for an expansion around the Penrose limit 
of AdS^ X S*^ in the context of the BMN correspondence |^), the above observations 
about the relation of these backgrounds to interesting space-time singularities provide 
additional impetus for understanding string theory in an expansion around such metrics. 

The analysis of the present paper can be generalised in various ways. In particular, even 
though we have only considered spherically symmetric metrics, the results we have ob- 
tained are certainly not restricted to spherical symmetry. A prototypical anisotropic ex- 
ample is the Kasner metric whose Penrose limit (together with those of other anisotropic 
or inhomogeneous cosmological models) was studied by Kunze ^H]- In many of these 
examples one sees that one once again gets a u~'^-behaviour either on the nose (for 
particular Penrose limits of the Kasner metric) or in the near-singularity limit. It would 
be interesting to establish a more general result along these lines, in particular also in 
view of the role played by the Kasner metrics (and Bianchi IX cosmologies) in the BKL 
discussion of the general solution of Einstein's equations near spacelike singularities pUj . 

Going beyond spherical symmetry may also shed light on a peculiar feature of the 
spherically symmetric metrics of power-law type studied in ^| ^] and here. Namely, 
we have seen that a rather special role is played by the "extremal" equations of state 
for which at least one of the inequalities in the DEC is saturated. As noted in this 
type of equation of state is not ruled out by physical considerations alone, and arises 
"rather too easily" near the spherically symmetric singularities of power law type that 
we have been considering. Indeed, as we have seen in section 4.4, whenever q > p + 2 
one obtains an extremal equation of state close to the singularity because then the 
"extremal" contribution to the curvature from the transverse sphere dominates the 
longitudinal (radial/time) contribution. Thus this may well be an artefact of spherical 

^This bound also plays an important role in a related context in the recent article \27\ analysing 
the propagation of scalar fields in space-time singularities, where the role of u is played by the usual 
tortoise-coordinate r, of Schwarzschild-like metrics. 
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symmetry. 

We have also not discussed null singularities in detail. To study these may not only 
be of interest in its own right but also because of the existence of a null counterpart 
of the BKL solution, a general solution of the Einstein equations describing a null 
weak singularity jHUI- There is now mounting evidence (see e.g. [311 132j ) that these 
singularities are realised inside realistic rotating black holes, where they are expected 
to form, as a consequence of the Penrose instability, at the location of what would have 
been the inner Cauchy horizon of the stationary Kerr metric. These null singularities 
bear some resemblance to the general rotating singular homogeneous plane waves found 
in jl5j . and one might wonder if Penrose limits of singularities inside rotating black 
holes give rise to these rotating singular homogeneous plane waves in the same way 
as rotating smooth (Godel) metrics give rise [231 the rotating smooth homogeneous 
plane waves of jl5j . For some comments on Penrose limits in this context see j34j . 

Since string theory is exactly solvable both in the non-rotating singular homoge- 
neous plane wave backgrounds J3] and for the general rotating smooth homogeneous 
plane wave metrics one also expects string theory in these general homogeneous 
singular backgrounds to be exactly solvable [HHI, and it will be interesting to find out if 
they actually arise from Penrose limits of realistic rotating black holes. 
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A Hamilton- Jacobi Equations, Geodesic Congruences and Adapted 
Coordinates 

Even though the Penrose limit itself only refers to a particular null geodesic, the original 
construction of Penrose limits in terms of adapted coordinates (Penrose coordinates) [HI 
|S] requires the embedding of this null geodesic into a twist-free null geodesic congruence. 

Indeed, in Penrose coordinates {U,V,Y^), i = l,...,d, the original space-time metric 
takes the form 

ds^ = 2dUdV + a{U, V, Y^)dV^ + 2hi{U, V, Y^)dVdY' + gij(U, V, Y^)dY'dY^ , (A.l) 
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corresponding to the {d + 2) coordinate conditions guu = gjjyk = 0, gjjy = 1. This 
exhibits the embedding of the original null geodesic V = = into a congruence 
of null geodesies parametrised by the coordinates V and Y^, the coordinate U being 
identified with the affine parameter of the null geodesies in the congruence. 

A systematic way of constructing Penrose coordinates makes use of the Hamilton- Jacobi 
(HJ) function of the null congruence jl2l I16j . In this Appendix, we shall explicitly 
construct the transformation from a general coordinate system to Penrose coordinates. 

The essence of the HJ method can be summarised by the observation that the momenta 

associated with the above null congruence ({/ = l,y = y^ = 0) are 

so that, in arbitrary coordinates x^, one has 

P^ = df,V . (A.4) 

Thus, since the geodesic congruence is null, g^^ d^VduV = 0, one can identify 

Vix^") = . (A.5) 

with the solution of the Hamilton-Jacobi equation 

5^^a^55,5 = . (A.6) 

corresponding to the null congruence 

x^' = g^'^dyS . (A.7) 

Conversely (see 1)2. 22() ) any solution S of the equations ()A.6IA.7|) gives rise to a (twist- 
free) null geodesic congruence and y = S is the corresponding null adapted coordinate. 

It thus only remains to understand how to construct the transverse coordinates Y^ . To 
that end we will now briefly review some facts about solutions to the Hamilton-Jacobi 
equations (see e.g. [3niEZl)- 

The general solution to the HJ equation HA.6|) can be rather involved but there usually 

exists a "complete" solution, complete in the sense that it depends on d + 2 integration 

constants |36| (d + 2 as in the rest of the paper is the space-time dimension).^ For 

^It is not always guaranteed that such a complete solution exists, though in all the cases that 
we consider here it does. The most general solution to (IA.6t is much more complicated and can be 
constructed from a complete solution by looking at x-dependent hypersurfaces in the space of integration 
constants by the method of envelopes |37| . 
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a complete solution to the HJ equation with integration constants a^, the associated 
geodesic congruence is = x'^(A, a^, Xq), where Xq are the positions of the geodesies 
at A = 0, Xq = x^(0, a^, j;g). The initial value surface parameterised by the Xq is a 
Cauchy surface for the HJ equation and can be represented algebraically by the equation 
F{x^) = 0. For a well-posed initial value problem, we require that the hypersurface 
F = has an everywhere timelike normal vector {dF)"^ < 0. 

One of the integration constants simply represents a constant shift of S. Furthermore, 
the HJ equation is homogeneous of degree two, so if S" is a solution, then kS, k = const ^ 
0, is also a solution. This scale invariance of the HJ equation is absorbed in the first 
order geodesic equations, ()A.7|) by a scale transformation of the affine parameter A. 
Therefore, there are only d non-trivial integration constants which we will denote by 
ak, k = 1, . . . ,d. 

Given a particular null geodesic 7, the integration constants Ok can be uniquely fixed. 
Indeed let p° = g^ux'^\x=o be the momentum of the geodesic 7 at A = 0. The mass-shell 
condition g^^p^Pu = is scale invariant and therefore there are d independent momenta. 
These can be used to determine the integration constants of the HJ function S via the 
equation 

pI = a^SU=o . (A.8) 

Therefore we can use the HJ equation to embed a given null geodesic into a twist free 
null geodesic congruence determined by the solution S. 

Given a null geodesic 7, the coordinate transformation from the original coordinates x^ 
of space-time to the Penrose coordinates can be defined using the HJ function S and 
coordinates Xq of the Cauchy hypersurface, as follows: 

We first parameterize the null geodesic congruence as described above 

x^ = x'^(A,x^) (A.9) 
F(x^) = . (A.IO) 

We have suppressed the integration constants because they are specified by the 
momentum of the null geodesic 7. Then we set 

U = \ 

V = 5K) . (A.ll) 

Note that S{x^) = S{xq) because 5" = 0. It remains to determine the coordinates 
yk {yovo. these data. For this observe that the level sets of S have null normal vector, 
because of HA.6|) . while the hypersurface F = Q has a timelike normal vector. Thus we 
have 

g^'''^^S^,F <^ , (A.12) 
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and the level sets of S intersect transversally the hypersurface -F = 0. The coordinates 
are found by solving the equations F{xq) = and S{xq) = V, i.e. Y'' are the 
coordinates of the transverse intersection of the -F = hypersurface with the level sets 
of the HJ function S. Using this, we can rewrite the first equation in ()A.10() as 

x^" = x^{U,x'^{V,Y'')) = x''{U,V,Y'') . (A.13) 

This is the transformation which relates a coordinate system on a space-time to the 
Penrose coordinates. 

Note that for a generic space-time, there is no natural choice for the hypersurface F = 0, 
i.e. for the function F. Instead F should be thought of as a gauge fixing condition which 
is chosen at our convenience. The Penrose limit metric does not depend on the choice 
of F, different choices simply corresponding to different ways of labelling the geodesies 
of the congruence on which the adapted coordinates are based. 

For the sake of completeness we will now show explicitly that in these coordinates the 
metric indeed takes the form HA.1|1 |12l I16j . First of all, we clearly have 

dx^ dx^ 

9UU = 9,.^^ = (A.14) 

because the geodesies x^(A,Xq) are null. Moreover, 

dx^^ dx" ^dx^^ dxf^ ^ ^ dV ^ /. -.^n 

guv = 9,.^^ = g.ug'-'dpS— = —d,S = ^ = l , (A.15) 

^^'^ _ dx^'dx" _ ^pdS dx" _ dV _ 

gui - g^^u^gyi - g^^ug a^^yi-^yi-O- (^-i^) 

B Generalisations: Brane Metrics, Isotropic Coordinates, Null 
Singularities 



It is straightforward to generalise the analysis of section 2.3 to include longitudinal 
world volume directions, 

f{r){-df)^ f{r){-dt^+dy^) . (B.l) 
A parallel frame in the brane worldvolume directions is Ei = f~^/'^dyi, and 

Bi,=dijdulogf{r{u))^/^ (B.2) 

which in turn leads to 

A^, = 5ijf{r{u))-y^dl\ogf{riu))y^ . (B.3) 
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The remaining of the components of A are as in section 2.3. 
Likewise, for isotropic coordinates, 

ds^ = -f{r)dt'^ + hir){dr^ +r^dnl) , (B.4) 

a straightforward calculation reveals that 

tr B = du log (^rr'^fh^ sin"'^^ (e) j (B.5) 

and 

^22 = 5„log(/i5rsin(0)) . (B.6) 

These lead to 

Bu{u) = du\og{rrhf^/^) 

Bd^) = d^f^dulog{rh'/^sm9) (B.7) 

with the corresponding second-derivative expressions for Aab{u). Again it is easy to 
include longitudinal directions. 

Finally, we consider spherically symmetric null metrics of the form 

ds^ = 2g{x)dxdy + /(x)dO^ . (B.8) 

The geodesic equations are 

x = Pg-^ , y = -^f-\ e = r'L, (B.9) 
where P and L are constants of motion. In this case, one finds 

tr B = du log{±gf sin(6')) (B.IO) 

and 

B22 = dulogip sm{9)) . (B.ll) 

Therefore, we have 

K2{u) = f{u)^/^sme{u) . (B.12) 
In particular, in terms of r(x) = /(x)^/^, A22 once again takes the standard form H2.40[) 

^22(^) = 44 - ^ • (B.13) 

r{u) r{u)^ 
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C Curvature of Szekeres-Iyer Metrics 



For reference purposes we give here the non-vanishing components of the Ricci and 
Einstein tensors of the metric (|4.5|) . 

ds'^ = rjxPdy^ - r^x^dx^ + x'^dflj (C.l) 

(for d = 2, these results can be inferred from Indices i,j refer to the metric gij of 

the transverse sphere (or some other transverse space) , with Rij and R the corresponding 
Ricci tensor and Ricci scalar. 

Ricci Tensor 

i(2p + 2qd + pqd - q^d)x~'^ 
\p{qd - 2)x~' 
R,j + \m{qd - 2)g,,x^-P-^ 

(d - + Irjqiqd - 2)gi,x'^-P-' (C.2) 
Ricci Scalar 

R = Rx-1 - ir/(4p + Aqd - d{d + l)^^)^-^^^^) 

= d((i - - ir/(4p + 4gd - d((i + (C.3) 

Einstein Tensor 

-\Rx-'i - \ridq[{d - l)q + 2p)x~^P+^^ 
-yid - - ^Tjdqiid - l)q + 2p)x-(P+2) 

-\Rx-'i + l'ndq{2p + 4 - (d + l)q)x-'^P+^^ 
-\d{d - + \rjdq{2p + 4 - (d + l)q)x-^P+^'^ 

G^x-'^ + hi^P - 4^ + 4gd - d{d - l)g2)5j2;-(P+2) 

-lid - l){d - 2)5ix-'i + \ri{Ap - 4g + Aqd - d{d - l)q^)6ix-^P^'^^ (C.4) 



Rx 



R, 



yy 



R. 



GZ 



Gl 



G] 
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